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1.  INTRODUCTION 


The  XBR-2D  code  is  •  two-dimensioiul  (2>D).  axisymmetric  beat  conductioii  code  qwciaUy 
fbnnulated  and  developed  Voitay  Technology  for  the  U^.  Aimy  Ballistic  Researeta  LjboiaiCHy  (now 
known  as  the  Anny  Research  Laboratoiy  [ARL])  to  simulate  gun  barrel  hett  tnnsfer  during  sing^  riiot 
and  burst  firing  scenarios*.  The  code  is  an  eiqdicit  fimte  difference  formulation  that  features  a  geometric 
radial  grid  netwoik.  This  provides  finely  placed  grids  at  the  bore  surface  to  accurately  model  foe  steq) 
temperature  gradients  at  foe  time  of  maximum  heat  flux,  while  defining  the  temperature  distribution 
forough  foe  barrel  wall  with  a  minimum  number  of  grids. 

The  oqdicit  tuture  of  the  formulation  requires  that  the  maximum  value  of  the  time  interval  be  limited 
by  qrecific  criteria  to  maintain  cmnputational  stability.  Thus,  the  maximum  time  step  is  omstrained  by 
foe  thickness  of  foe  smallest  grid. 

The  burst  fire  event  amsists  of  a  number  of  brief  ballistic  cycle  periods,  vriiidi  are  generally  less  foan 
1(X)  ms  in  duration  and  are  diaracterized  by  intense  heating  and  steq>  temperature  gradients  at  the  bore 
surface.  These  represent  foe  most  severe  amditions  that  foe  code  is  designed  to  simulate.  Between  foese 
periods  of  ballistic  activity  are  dwdU  periods  in  which  beat  transfer  is  low  and  foe  tbrou^-wall  thermal 
profile  relaxes  through  heat  conductkm  processes.  The  temperature  gradient  at  the  bore  surface  decreases 
in  magnitude  but  the  grid  pattern  remains  extremely  fine.  Prior  to  this  effort,  foe  same  time  step  and  grid 
pattern  were  used  foroughout  foe  simulatitm  to  evaluate  both  the  brief  ballistic  cycle  and  foe  potentially 
extoisive  dwell  period  between  shots.  When  simulating  long  burst  fire  sceruuios  for  artillery  with  the 
original  version  of  foe  XBR>2D  code,  foe  computer  time  is  excessive;  execution  times  of  20  hours  on  foe 
ARL’s  Cray  computer  have  been  cited.  Thus,  it  is  deared  to  reduce  the  run  time  exhibited  by  foe  code 
while  maintaining  its  computational  accuracy. 

1.1  Approaches  to  the  Problem.  Two  possible  approadies  were  evaluated  to  reduce  foe  executkm 
time.  The  first  rqrproach  involved  changing  the  numerical  method  to  either  an  implicit  finite  differetKe 
or  a  boundary  element  method.  An  implicit  finite  difference  scheme  involves  solving  foe  entire  grid 
matrix  simultaneously  for  each  time  step  as  opposed  to  foe  explicit  scheme  which  solves  each  grid 
indqxndently  based  on  the  values  at  foe  previous  time  step.  An  advantage  of  foe  implicit  method  is  that 
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there  are  no  stability-related  time  step  reactions.  The  main  constiaim  on  the  time  step  is  the  need  to 
maintain  adequate  resolution  of  the  heat  trensler  time  histoiy.  During  the  dwell  period,  the  step  can  be 
increased  without  sacrificing  accuracy.  Using  test  cases  designed  to  evaluate  the  difierent  mohods,  die 
im{dicit  finite  difference  solutions  did  not  agree  with  the  ex{dicit  method  in  any  calculatkms.  Afim 
reviewing  both  methods  to  verify  that  diey  were  {Mopeily  implemented,  and  calculating  analytical  solutions 
to  die  test  cases,  it  was  detennined  diat  the  imi^cit  method  was  not  providing  accurate  solutums  to  the 
proUems.  so  it  was  eliminated  as  an  opdtm. 

The  boundaiy  element  method  (BEM)  has  the  possibility  of  significandy  improving  the  qieed  as  well 
as  the  accuracy  of  the  sdudoa  The  BEM  solves  the  governing  equadmi  cmly  along  the  boundaiy. 
eliminating  the  need  to  compute  internal  points.  If  internal  values  are  desired,  rady  the  poim  of  interest 
needs  to  be  calculated  This  eliminates  the  large  number  of  extraneous  internal  points  which  need  to  be 
calculated  in  either  of  the  finite  differoice  methods.  However,  the  BEM.  e^iecially  for  the  case  of 
tiansiem  heat  conducdon.  involves  very  complicated  mathematics  in  its  derivation.  Because  of  the  limited 
scope  of  this  task  order,  the  BEM  could  not  be  pursued  to  its  conclusion,  although  smne  very  valuable 
and  promising  woik  was  accomidished. 

The  second  a|iproach  to  speeding  up  the  XBR-2D  code,  and  die  one  that  was  ultimately  implemented, 
was  to  woik  with  the  current  eiqilicit  finite  difference  scheme  and  improve  its  efficiency.  Since  the  time 
stq)  in  the  original  versimi  of  die  code  was  constant  for  both  die  ballistic  cycle  and  the  dwell  period,  an 
increase  in  the  stq>  size  for  die  dwell  period  (when  the  stability  requirements  are  rHaxed)  allows  for  a 
decrease  in  the  execution  time.  An  algorithm  was  developed  which  computes  die  maximum  temperature 
gradient  in  the  system,  and  dynamically  alters  the  time  st^  based  on  this  value.  This  algorithm,  alcmg 
with  other  programming  optimizations,  decreased  the  execution  time  by  factors  of  diree  to  six,  depending 
upon  the  loigth  of  die  dwell  period.  These  improvements  in  the  execution  time  were  achieved  without 
significantly  affecting  the  numerical  results  of  the  XBR-2D  simulations. 

2.  EVALUATION  OF  NUMERICAL  ALTERNATIVES 

To  quickly  evaluate  the  various  numerical  methods,  two  test  cases  were  devdoped.  Case  one, 
illustrated  in  Hgure  1.  is  a  fiat,  two-dimensional,  6x6  plate  with  a  constant  temperature  of  300°  on  one 
side  and  a  constant  temperature  of  1.500°  rai  the  o^xisite  side.  The  remaining  sides  are  insulated.  The 
steady-state  solution  to  diis  problem  is  a  linear  distribution  fiom  300  to  1.500. 
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Hie  second  leA  case  is  also  Ulustnted  in  Hguie  1.  In  diis  case,  die  plate  fiom  die  fiist  test  case  is 
insulated  on  three  sides  and  a  heat  flux  is  applied  to  the  fiiuith  side.  The  heat  flux  is  qiplied  for  30  ms 
and  is  followed  by  a  pause  of  60  ms.  The  cycle  is  then  repeated  to  roughly  simulate  die  oonditions  in  a 
gun  barrel. 
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Figure  1.  Heat  transfer  test  cases  used  to  evaluate  various  numerical  tedmioues. 


2.1  Explicit  and  Implicit  Finite  Difference  Methods.  The  numerical  method  contained  in  XBR-2D 
version  1.0  is  an  explicit,  centered  finite  difierence  method.  This  method  is  explained  in  detail  in  the 
original  XBR-2D  report  (Chandra  1990).  In  an  x-y  coordinate  system,  the  transient  heat  equation. 


A*T 


I  6T 

oTT  ’ 


(I) 


can  be  modeled  using  this  method  with  the  centered  difference  equation 


^  1  ««  B 


*  oAt 


Tuu  -  2T“,  *  T“,.,  ^  T,V,  -  2T,:i  *  T.^j., 


Ax* 


Ay* 


+  At^ 


(2) 


AX 

0  >  -k  ^  s  Boundary  Condition 
Ax 


n^re  ^  is  the  external  heat  flux  at  the  grid,  Ty  is  the  temperature  of  the  grid  of  interest,  and  Ty+I,  Ty.j, 
Tmj.  T|.tj  ate  die  temperatures  of  adjacent  grids. 
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The  tiiwg  ttep  was  varied  with  die  exfdidt  method  until  a  converged  solution  was  reached.  This 
process  is  fllustiated  in  Figures  2~4.  Figure  2  shows  the  overall  behavior  of  die  first  test  case  widi 
various-sized  time  steps.  Temperature  differences  widi  the  larger  time  steps  all  occur  in  the  very  eariy 
r^kn.  Hgures  3  and  4  fiicus  tm  die  early  time  behavior  where  the  solutions  cmiveige  as  the  step  siz? 
is  reduced.  A  time  stqi  of  0.010  s  was  diosen  for  die  baseline  calculations. 

The  results  of  die  second  test  case  are  shown  in  Rgure  S.  For  this  case  the  exfdidt  method  was 
unstable  if  the  time  stqi  was  large,  but,  (mce  die  stqi  was  small  enou^  to  maintain  stability,  die  results 
remained  the  sane  as  the  time  step  was  reduced  further.  A  step  size  of  0.001  s  was  chosen  for  the 
baseline  calculadons  of  Test  Case  No.  2. 

These  baseline  solutions  were  then  used  to  evaluate  the  accuracy  of  die  implidt  finite  difference 
calculations. 

The  imidicit  finite  difference  method  chosen  to  be  evaluated  was  die  alternating-direction  imjdidt 
(ADI)  method  (Anderson,  Tannehill,  and  Pletdier  1984).  This  method  was  chosen  over  standard  centered 
difference  implidt  methods  because  it  produces  a  tridiagonal  system  of  equations.  A  tridiagmial  matrix 
requires  substantially  less  computational  time  to  solve  dian  nontridiagonal  matrices  sudi  as  diose 
developed  with  standard  implidt  methods.  The  key  to  minimizing  die  execution  time  of  die  implidt 
method  is  die  sdection  of  an  algorithm  to  solve  the  tridiagtmal  system  of  etpiations.  Three  tridiagnial 
matrix  solving  algorithms  were  evaluated:  the  Thcxnas  algorithm  (Anderson,  Tannehill,  and  netcher 
1984),  Gauss-Jordan  Elimination,  and  LU  Decompc^itimi  (Press  et  al.  1989).  Hie  Thranas  algoridim, 
found  to  be  slightly  faster  than  die  other  two,  was  chosen  as  the  solution  mediod.  All  three  algorithms 
gave  identical  numerical  results. 

Figures  6  and  7  show  a  cmnparistm  between  the  erqilidt  and  imididt  solutions  for  die  two  test  cases. 
The  figures  show  a  distinct  difference  in  the  behavior  of  the  implidt  method.  Both  finite  difference 
methods  were  re-examined  and  rqspeared  to  be  implemented  properly,  but  the  solutions  could  not  be  made 
to  agree.  An  analytical  solution  was  developed  for  the  second  test  case.  When  diis  solution  was 
cmnpared  to  die  implidt  and  exididt  results,  it  became  obvious  that  die  imididt  method  was  not  {ntqieriy 
solving  the  proldem.  The  irrqdidt  method  was,  therefore,  diminated  as  an  optiort 
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Figure  4.  Eariv  time  behavior  for  various  in  the  exoBcit  Test  Ctee  No.  1. 
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2^  Boundafv  Bemew  Method  (BEM).  As  the  name  impUes,  in  the  BEM  approach  the  govemmg 
dilfeieotial  equatioiis  are  trensfoimed  into  integral  identities  \KMdt  are  qtfdied  on  the  boundary.  These 
initials  are  dien  numerically  integrated  over  the  boimdary.  whidi  has  been  divided  into  small  segmeitts 
(boundary  deoMnts).  The  BEM  cm  be  applied  to  comidex  boundary  shares  relatively  easily  and. 
ftinfaennore.  since  all  die  approximations  are  restricted  to  the  boundary,  it  cm  model  regioos  having 
r^ridly  chmging  variables  with  considetable  accuracy. 

Since  the  problem  is  solved  only  on  the  boundary,  and  fewer  calculations  are  involved,  die  BEM 
should  provide  a  significant  decrease  in  the  execution  time  ova*  eidier  finite  difference  mediod.  Fbr  die 
Poisson  eqtuoion,  Q.e.,  steady  state  heat  omduction),  the  BEM  formulatim  is  well  known  and  documented 
(Becker  1992;  Btebbia  and  Dominguez  1989;  Biebbia  and  Walker  1980).  However,  when  a  time 
dependency  is  added  to  the  (»oUem,  the  formulation  beomnes  much  more  complex. 

An  existirig  steady  state  heat  omduction  BEM  code,  POLINBE  (Brebiria  and  Dmninguez  1989),  was 
chosen  fix*  evaluatioiL  This  code  was  written  to  model  proUems  with  geometries  similar  to  die  test  cases 
used  to  evaluate  the  finite  difference  techniques.  When  run  as  a  steady  state  proUem.  a  solution  was 
obtained  practically  instantaneously. 

Appendix  A  describes  in  detail  the  mathematical  formulation  of  die  transient  heat  conduction  BEM. 
This  formulation  was  added  to  die  POLINBE  code  in  m  attempt  to  develqi  a  transient  BEM  solution  to 
model  die  test  cases.  Several  proUems  where  encountered  in  the  formulation  of  the  mediod.  The  most 
important  difficulty  was  the  selection  of  m  qipioximatitHi  functioa  This  functim  is  a  key  demem  of  die 
tranaem  heat  conduction  foimulatim  and  further  effort  is  required  to  determine  an  iqipropriate  function 
for  the  band  heating  problem. 

The  BEM  is  quite  complex,  and  devdopmem  and  implementatim  of  this  teduiique  could  not  be 
acemnplished  during  die  presort  program;  however,  the  progress  that  was  made  strongly  suggests  that  it 
holds  great  potential.  It  is  bdieved  that  this  method  may  be  able  to  provide  a  substantial  decrease  in 
execution  time,  while  providing  an  increase  in  computational  accuracy;  however,  these  attributes  have  yet 
to  be  demmstrated  for  the  band  heat  conduction  problem. 
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3.  XBR-2D  CX)DE  MODinCATIONS 

3.1  Time  Step  Changes.  The  test  cases  described  in  Section  2  iUustnted  the  difficulties  asaodMed 
widi  adiieving  execution  time  improvements  by  altering  the  integration  nmthod.  The  sample  cases  were 
not  able  to  maintain  a  Consistent  numerical  stdution  without  additional  effim  that  would  have  lequiied 
giemeriesouroes  than  were  available.  It  was  decided  that  the  best  option  m  this  time  was  to  improve  the 
existing  explicit  finite  difierence  method’s  efficiency  modifying  the  time  step  and  optimizing  die 
FORTRAN  coding. 

In  die  ori^nal  version  of  die  XBR-2D  code,  die  time  step  (At)  was  based  on  the  stability  criteria: 


At  « 


Ar* 

6a 


(3) 


Fbr  a  barrel  with  a  duome  {dating.  Ar  was  fixed  at  the  {dating  thickness,  otherwise  Ar  was  determined 
based  on  a  user-{HDvided  At  This  time  step  is  dim  used  throughout  the  ballisdc  cycle  and.  in  die  case 
of  mult^e  shots,  die  dwell  period  between  shots.  In  cases  where  there  is  a  significant  dwell  period,  die 
tem{)erature  gradients  can  becmne  quite  small  and  a  step  as  small  as  dictated  by  equation  3  is  not 
necessary  to  maintain  stability. 

In  a  first  attempt  to  ^reed  tq)  the  code,  die  time  sxep  was  increased  equaUy  for  the  endre  dwell  period. 
Any  significant  increase  in  die  At  caused  the  code  to  become  unstaUe.  To  maintain  stability,  the  time  step 
modification  was  then  changed  to  a  function  depeadeOL  o{X>n  the  maximum  radial  tem{)er8ture  gradient 
in  die  system.  In  the  dwell  region,  the  stq>  is  defined  as 

At  -  CF  (4) 


CF  is  a  geometric  ctxiection  factor  whidi  accounts  for  differenoes  in  grid  sizes  between  various 
gun  calibers.  At^  is  defined  in  equatitm  3.  and  (VT,)y.  is  the  maximum  temperature  gradient  in  die 
band.  The  dwell  time  stq>  is  restricted  so  that  it  cannot  be  smaller  than  the  base  time  stqi.  At^  nor  can 
it  be  greater  than  1(X)  times  die  base  time  stqi. 
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This  condlsrion  picfvkles  die  Ivgest  increase  in  step  tize  when  die  heat  conduction  approadies  a 
steady  stale,  and  will  provide  the  most  sidistMttial  speed  increases  in  multiple  dmt  simulatians  where  die 
kngdi  of  die  dwdl  peiiod  is  long  oonqiared  to  the  length  of  the  baUisdc  cyde,  as  in  die  case  of  aitiUety. 

3^  Code  Optimization.  Many  small  code  optimizations  were  also  made  in  most  of  XBR-2D’s 
subroutines.  These  involved  iqdachtgfiequendy  used  anay  elements  (sudi  as  T(l't-l))widi  local  variables 
(audi  as  Tl),  iqdadng  vqiedtive  calculations  with  a  single  calcdadon,  and  moving  oeitain  calculatkms 
oittofkMps.  These  opriisizatioos  were  most  significam  in  the  subroutine  HEATRANSudddicairies  out 
die  finite  difference  cateuladons.  and  in  the  subroutines  which  calculate  the  boundaiy  conditions.  These 
modifications,  on  their  own,  decreased  die  execution  time  of  XBR-2D  by  15%. 

An  option  diat  was  induded  in  the  original  XBR-2D  was  the  ability  to  turn  off  the  screen  iqxlates. 
By  mining  off  die  screen  iqxlates  alone,  the  execution  time  on  a  PC  decreases  by  15%.  For  die 
comparisons  done  in  the  next  section,  die  screen  update  was  mined  off  for  both  the  original  XBR-2D  tuns 
and  the  new  runs.  All  of  die  runs  in  die  next  section  were  calculated  on  a  33~MHz  80486  PC-compatible 
con^Mter. 

4.  CODE  EVALUATION 

4.1  155-mm  Howitzer.  A  155-mm  howitzer  test  case  for  Unichaige  was  provided  by  the  ARL.  This 
case  simulates  die  firing  of  60  rounds  with  a  5-s  dwell  time  between  each  shot  Initially,  to  qieed  iqi  the 
debugging  and  testing  process,  only  die  first  round  was  studied.  This  case  involved  a  single  shot  followed 
by  a  dwell  of  5  s.  Table  1  coaq»fes  die  bore  suifoce  tempenture  rise  of  diis  single  dwt  case  before  and 
after  the  code  modifications.  TtMe  2  compares  the  peak  bore  suifroe  temperature  at  each  locatioa  The 
values  are  unchanged  for  the  two  simulations.  The  modifications  to  the  code  do  not  take  effect  umil  after 
the  peak  temperature  has  been  reached.  The  new  code  required  less  than  20%  of  die  mtecution  time 
required  by  veiskm  1.0  of  XBR-2D. 

Figure  8  is  a  grqih  of  the  probe  iemperaturc  for  tiie  new  and  old  nins.  The  two  cases  are  so  dose 
dutt  any  differences  are  riot  discemitale  on  die  graph. 
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Tride  1.  Cba^Mriaon  of  Old  and  NdvXKt-2D  Suffice  Teaipeniuie  Rite  fx  a  Sio^Shot 
155-auB  Teat  Caae  (Ten^maittiea  are  at  the  Gompireion  of  tte  Simulation) 


Run 

Eaecution 

Time 

(a) 

Temperature  rise  (K)  at  each  location. 

Locations  in  inches  finm  die  breedi  end. 

lao 

42.3 

85.7 

129.1 

172.5 

215.9 

259.3 

310.8 

Oiig 

3.6S0 

29.90 

27.91 

22.32 

20.73 

18.70 

15.74 

16.81 

15.81 

New 

705 

29.75 

27.91 

22.32 

20.72 

18.70 

15.74 

16.81 

15.81 

Table  2.  Cbmpaiiaon  of  Old  and  New  XBR-2D  Peak  Surface  Temperature  Rise  for  a  Single-Shot 
ISS-mm  Test  Case 


Run 

Execution 

Htne 

(s) 

Peak  tenfietature  rise  (K)  at  each  locatioa 

Locations  in  inches  firom  die  breech  end. 

10.0 

42.3 

85.7 

129.1 

172.5 

215.9 

259.3 

310.8 

Otig 

3.650 

883.3 

1031.7 

847.3 

695.0 

563.8 

418.1 

443.9 

357.3 

New 

705 

883.3 

1031.7 

847.3 

695.0 

563.8 

418.1 

443.9 

357.3 

Table  3  summarizes  the  results  of  die  foil  60-iound  simulation.  Hguies  9  and  10  are  plots  of  die 
surface  tempoature  rise  and  the  probe  temperature  rise,  reflectively.  There  are  no  discemitde  differences 
ondie  grqihs. 


Table  3.  Comparison  of  Old  and  New  XBR-2D  Surface  Temperature  Rise  for  a  60-Shot  ISS-mm  Test 
Case  (Temperatures  ate  at  die  Comptetion  of  die  Simulation) 


Execution 

Hme 

(8) 

Temperature  rise  (K)  at  eadi  locatioa 

Locations  in  indies  from  the  Ixeedi  end. 

10.0 

42.3 

85.7 

129.1 

1723 

215.9 

259.3 

310.8 

Orig 

264319 

3733 

368.0 

3122 

275.4 

242.3 

202.4 

213.7 

202.0 

New 

78,144 

3723 

367.7 

3113 

274.7 

241.8 

202.0 

213.3 

201.6 
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Hgure  8. 


140 

120 

100 

60 

60 

40 

20 

0 


Figiiie9. 
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4^  iM-irnn  Tank  Gun.  A  test  case  modeling  a  120-mm  tank  gun  was  assembled  by  Veritay.  A 
duee-sbot  bunt  of  M829  tank  rounds  was  modeled  with  a  10-s  dwell  time.  The  peak  smftoetenqwntuies 
are  summarized  in  Table  4.  Figure  11  shows  die  surface  temperature  rise  for  die  diird  shot  of  die  burst 
The  new  values  were  obtabned  after  die  code  coirectkms  described  in  Appendix  C  were  implemeDied. 

4.3  7<t4nin  riaion.  A  test  case  modeling  a  76-mmOtoMelara  gun  was  assemUed  by  Veritay.  This 
shnularion,  summarized  in  Table  S,  included  a  duee-round  burst  firing  with  a  1-s  dwdl  time  between 
diots.  hi  this  case,  die  new  XBR-2D  stanulatkm  was  four  times  faster  duni  the  original  XBR-2D 
rimoladon.  Figure  12  compares  die  suifoce  temperature  Mstoiy  at  die  22-in  location  and  includes  a 
magnified  view  of  the  diird  temperature  peak.  The  pointed  shape  in  the  "old-XBR"  curve  is  a  result  of 
plotting  only  every  tenth  point  The  new  values  were  obtained  after  the  code  correcdons  described  in 
^jpendix  C  were  imjtonented. 
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Table  4.  Comparison  of  Old  and  New  XBR-2D  Suiface  Temperature  Rise  for  a  Single-Shot 
120mm  Test  Case 


Run 

Executkn 

Tune 

(8) 

Suriare  ieaqxmure  rise  (K)  at  each  location. 

Locations  in  inches  from  the  breech  end. 

13 

20 

23 

33 

40 

43 

30 

Orig. 

15^39 

1,032S 

1,119.1 

1,206S 

1^.6 

1330.7 

13124 

13033 

1394J0 

New 

1^647 

1^1.7 

1.119J0 

1,208.1 

1^408.6 

13384 

13233 

13213 

1312.7 

Figure  11.  Surface  temperatures  at  the  IS-in  location  for  flic  120-inm  gun  simulations. 
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Tabk  5.  Oompafison  of  Old  and  New  XBR*2D  Suiftoe  Tempoaiure  Rise  Ah’  a  Single-Sbot 
76-iniii  Test  Case 


Run 

Execution 

Time 

(8) 

Peak  surface  temperature  rise  (K)  at  eacli  locadoa 

Locatimis  in  indies  from  die  breedi  end. 

22.0 

30.0 

40.0 

50.0 

60.0 

70.0 

Oiig. 

4.4S2 

883.4 

86S.S 

746.5 

611.4 

510.3 

539.8 

New 

1.1S8 

864.9 

871.2 

762.6 

631.9 

5252 

560.4 

0  500  1000  1500  2000  2500  3000 

Time  (ms) 


I^guie  12.  Surface  temperature  history  at  the  22-in  location  for  the  76-mm  am  simulation. 
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5.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  XBR-2D  code  was  modified  and  execution  time  was  decreased  by  a  factw  of  52  in  the  IS5-mm 
howitzer  sample  case  provided  by  the  ARL.  The  decrease  in  execution  time  was  mostly  obtained  by 
modulating  the  time  stq)  during  the  dwell  period  between  shots.  The  Itmger  the  dwell  period  relative  to 
the  ballistic  cycle,  the  greater  the  speed  imiuovemait  will  be.  Several  programming  optimizatians 
provided  additional  increases  in  efficiency  and  speed. 

Techniques  for  adapting  BEMs  to  gun  barrel  heating  inoblems  were  explored.  The  method  was  used 
far  a  steady-state  heat  omiductitm  problem  and  it  tqrpears  very  {nomising  for  both  its  accuracy  and  its 
speed.  The  formulation  of  a  transient  heat  transfer  BEM  is  a  very  complex  process,  and  the  limited  scqpe 
of  this  effort  did  not  allow  this  work  to  be  com|deted.  The  theoretical  background  was  developed  to 
implement  this  mathematical  technique  in  solving  the  2-D,  axisymmetric  barrel  heating  problem;  however, 
several  obstacles  remain  that  prevented  its  adoptitm  in  diis  program. 

The  Army  may  wish  to  consider  funding  an  additional  effort  to  incorporate  this  state-of-the-art 
numerical  method.  This  effort  would  include  development  of  approximation  fiinctiotrs  required  to 
optimize  the  BEM  technique  for  the  transient  barrel  heating  problem,  and  permit  evaluation  of  the 
technique  in  terms  of  computational  accuracy  and  speed.  Another  item  that  would  improve  the  accuracy 
of  the  computed  results  is  the  indusitm  of  temperature-dependent  material  properties.  At  presort,  thermal 
diffesivity  and  thermal  conductivity  are  treated  as  constant  values  in  the  formulation  and  integration  of 
the  governing  equations. 
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APPENDIX  A: 

BOUNDARY  ELEMENT  METF^  J 
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BACKGROUND 


The  boundaiy  donett  metliod  (BEM)  has  received  cooiiderdite  attettion  over  tlie  past  15  yem.  and 
qjpean  to  hold  mudi  psomiae  for  use  in  adving  transieot  heat  conduction  profaleiDS,  alttiou^  ceitdn 
dilBculttes  renudn  to  be  adved. 


The  BEM  ^ipareidy  originated  as  a  adieme  for  numerical  computations  laigely  as  an  extension  of 
approximation  tedmiques  invdving  Green’s  functions,  int^ral  equations,  finite  element  sdiemes,  and 
finite  difference  qtptoaches  (Bidibia  and  Walker  1980).  The  initial  qtplications  of  BEM  were  directed 
toward  time-indqiendent  qtatial  problems,  sudi  as  steady-state  heat  fiow,  potential,  and  elastostatic 
problems.  The  apfdicatkm  to  time-dependent  fields,  including  transient  heating,  has  required  additional 
attentitm  to  detennine  effective  means  of  dealing  with  the  time-dqxndenoe  feature.  This  activity  is  still 
on-going. 

A  numbo-  of  different  BEM  fonnulati(ns  have  been  advanced  for  treating  transient  heat  conduction 
problems.  The  most  successful  from  ttie  stanfooint  of  accuracy  is.  according  to  Wrobd  (1988),  the  one 
vriiidi  enqdoys  time-dqtendent  fundamental  solutums*.  This  approadi  is  a  numoical  extension  of  foe 
classical  Green’s  function  mefood  (Morse  and  Feshbadh  1953;  CZanflaw  and  Jaeger  1948)  for  obtainiog 
analytical  stdutions  to  simple.  tractaUe  difiusion  inoUems.  Since  no  mefood  has  yet  been  found  to 
conqttetdy  sqMiate  foe  time  dependence  expUddy  fiom  foe  unknowns  in  foe  BEM  equations  (even  when 
foe  time-dqiendent  fundamental  solution  is  used),  it  is  still  necessary  to  use  a  time  marcbing  sdieme  of 
some  type  to  obtain  a  numerical  solution  of  foe  BEM  equations. 

Two  such  time  marching  schemes  have  typically  been  used.  In  the  first,  foe  total  time  is  divided  into 
stqts  with  eadi  new  stq)  treated  as  a  new  problem.  Wifoin  foe  initial  time  stqp  foe  initial  and  boundary 
conditions  can  be  used  to  solve  foe  problem  for  foe  unknown  variaUes.  The  temperatures  at  a  number 
of  ptrints  wifoin  foe  domain  are  calculated  and  used  as  initial  values  for  foe  second  time  stq).  together  wifo 
foe  origitud  boundary  values,  to  ^ve  the  problem  anew  for  the  secrnid  time  interval  This  process  is 
repeated  stq>  by  stq>  until  the  final  time  is  reached.  The  solution  may  become  unstaUe  if  foe  time  stq)s 
are  very  small  (Becker  1992). 


IliB  tamiBolDgy  "findanMital  lolutian,"  nforrins  to  dto  govemais  pvtul  differantial  equiiiHW,  it  a  tohiikm  ofatanied  for 
die  infiiMie  AMMin  widtoot  ttldQg  into  conakfaruioii  qiecific  finite  bounday  conditiont  of  die  probfem. 


21 


In  dK  Moond  miiching  scheme,  the  int^ratiaos  widi  re^Kct  to  time  always  stan  from  the  iniiial 
conditions  and  die  effect  of  die  time  evoludon  of  heat  tiansfened  is  detennined  duou^  boundaiy 
im^ndons,  so  that  values  of  teoqierBture  at  ptdms  widdn  the  doouin  of  interest  are  not  lequiied. 
However,  the  atdudon  at  ny  dme  stq>  is  sdn  dqxndent  on  the  stdudon  at  all  previous  time  sfiqie  (not  just 
the  previous  one  as  in  die  first  marching  sdieme)  going  back  to  the  initial  dme.  This  process  suffers  ftom 
requirements  for  huge  amounts  of  computadon  dme  and  storage  tpKt  (Becker  1992). 

From  a  computational  stam^iomt.  these  approaches  using  the  time-dependent  fundamental  solution  are 
reoogniaed  as  having  significant  limitadons  (Wrobel  1988). 

Ahemaie  BEM  fonnuladons  of  die  transient  heat  conducdon  protdems  include  the  use  of  Ijqilaoe 
tnnsfonns  to  remove  the  dme  dependency,  solve  the  equations,  and  then  apply  die  inverse  Laplace 
tiansfonn  to  obtain  physical  values  (Brehbia  and  Walker  1980;  Rizzo  and  Stuppy  1990).  Another  involves 
a  coiqded  boundaiy  element-finite  differmce  (BE-FD)  iptnoadi  in  udiidi  finite  difference  dme  steps  are 
used  (CUiran,  Lewis,  and  Gross  1980).  This  coufded  sdieme  is  similar  to  the  first  dme  marching  sdieme 
noted  above.  A  third  approadi  uses  a  potential-type  method.  This  is  an  indirect  version  of  the  time- 
dependent  Green’s  foncdon  mediod  (Cunan,  Lewis,  and  Cross  1986).  Several  odier  approadies  are  also 
known  (Wrobel  1988). 

In  recent  years,  diere  has  been  a  trend  towards  the  use  of  time-indqxndent  fundamental  solutions 
togedier  with  an  approximatB  treatment  of  die  time  derivadve  term.  Two  key  references  for  this  approach 
include  die  wodts  of  Wrobel  Brebtna,  and  Nardini  (1986),  and  Wrobel.  Telles  and  Brdibia  (1986). 
Claims  of  significant  reductions  in  computadon  time  and  storage  requirements  have  been  made  vdien  usmg 
this  qiproach  (Wrobel  1988;  Wiobd,  Telles,  and  Brebbia  1986). 

Although  a  numbw  of  odier  BEM  foimuladons  were  examined,  efforts  here  enqdiasized  dds  time- 
indqiendent  fundamental  solution  ^iproach.  The  mathematical  formulation  of  diis  BEM  qiproadi  is 
sununarized  in  the  folfowing  section. 


22 


tlieaNDEffiNDem' FUNIMMEOTAL  SOUmON  BEII 


lliB  tnmkat  hBtt  cquirion  to  be  srived  for  the  temperamre  dtotribudon  ii(M)  over  ite  domiin  Q  if 


V*u  (x.t)  -  I  4  (*.t)  . 

K  M 


(A-1) 


wbeie  x  is  a  fidd  point*  with  oooidinaies  (x,.  x,,  x^  in  Q  or  on  its  boundaiy  F.  and  k  is  a  (oonstant) 
difliisivity  coeffidenL 

Boundaiy  oonditiaQs  aie  taken  heie  to  be  of  die  simple  fimn 


u<x)  -  u(x)  on  r,  , 


(A-2) 


wtaeie  r  «  r,  -t-  r,  is  die  complete  boundaiy  of  Q,  q(x)  is  the  flux,  n(x)  is  die  unit  outwaid  nonnal  at 
point  X.  md  k  is  the  (constant)  diennal  oonducdvity  coeSicienL  These  boundary  conditions  have  been 
extended  to  indude  nKfladan  effects  by  Dd^gueiiedo  and  BidMa  (1988)  and  Onialii  and  Knioki  (1982). 

Ihe  initial  condition  at  time  t  s  o  is  of  tiie  form 

u(x.o)«tt,(x)  on  .  (A-3) 

Now  apply  Gieen’s  second  identiqr  to  a  pair  of  die  functions  u  and  u*.  udieie  u*  is  die  fundamental 
solution  to  La(daoe*s  equation  (not  die  tiansient  heat  equation) 

kV>u*(C.x)- -S‘(C.»)  .  (A-4) 


'nw  lann  "Said  point*'  in  die  ament  BEM  Itewin  lefoe  lo  •  leaeral  nriefale  point,  edmeoc  "aauoe  pomt”  (lo  be 
i**"^**^  in  die  foUiiwing)  ie  tehm  to  be  e  find  point  for  eny  one  celcaletion. 
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and  i^ieie  8^  ftjO  te  Dine  ddta  Itaiictk^  The  ftmctkio  u  iqiitaentt  the  Sdd  genenied  by  a 
oonccniniBd  unit  aomce  acting  at  point  4.  lUt  ^vea  as  a  result 


J(u*kV*u  -  ukV*u*)  dQ  ■  j|(u*q  -  q*u)  dT 


(A.5) 


udiere 


(A-d) 


Recall  that  in  equation  A-S,  u  and  q  are  ftnetioos  of  x.  u*  and  q*  are  Ihnctions  of  ^  and  x,  and  the 
integntions  are  canied  oat  with  reqrect  to  the  vailtibte  x. 

Subatitoting  equation  A*4  into  A-5  gives  the  result 


^TtV^idD  ■  -  u' 


jffu'q  - 


q*a)  dT 


(A-T) 


irfme  u'  is  the  value  of  u  re  the  aouioe  point  For  my  other  fixed  position  of  die  unit  aouice  i  which 
might  be  chosen,  the  coneyonding  values  of  u*  and  q*  will  be  different,  and  a  new  integml  equation  (of 
identical  fonn)  will  apply. 

The  fundamental  setiution  of  equation  A-4  fiv  an  isotropic  three-dimenaioaai  (3-D)  domain  is 

u*-^  .  (A-8) 

4«r 

and  fire  a  two-dimensional  (2-D)  domain  is 

whree  r  is  die  fiom  die  unit  point  aouioe  re  ^  to  ai^  odier  point  under  consideration. 
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Bqoaiaa  A<7  is  vaM  for  wy  find  poinc  4  in  ibe  dooutau  bnt  10  fooNitaie  the  pfobicm  ai  a  boundiiy 
ieGtaiqie,tfiB|ioitt(anitt  teuton  10  the  boonduy.  Upon  axouating  for  ifaif.  wi  for  the  juny  <rf  the 
taiepil  in  <{*.  fliB  IbOowinf  Inieinl  eqn«kn  is  otoined  (Wiobd.  Telks,  and  BvdMa  1986): 


•  -  c‘a‘  ♦  J  (tt*q  -  q*tt)  dT  (A-IQ) 

where  c*  is  a  ftmction  (tf  the  solid  aigie  the  bouodaiy  makes  at  a  point  ^  (see  Brebbia  and  Dnmingues 
1989,  p.  7Q),  sod  the  taiegnl  in  q*  is  calCBlaied  in  die  Caudiy  princ^de  vaine  sense. 

In  the  2-D  case,  Cor  example,  if  die  boudaiy  is  smooth  at  the  point  4  diere  is  no  disconthwigr 
ofdope),  the  value  of  c*  is  M  father  than  unity  when  4  is  an  intBfior  point  If  a  disooodnuity  of  slope 
occnis  at  4  dien  cl «  8/2  x,  where  6  is  the  imenial  angle  of  the  boundary  F  at  The  quantity  c*  applies 
to  the  continuous  boundaiy.  as  well  as  to  die  points  udme  dements  of  die  boundaiy  join  when  the 
boundaiy  is  divided  into  a  finhe  nnmber  of  dements  for  numeiicd  andysis. 


PUither,  substituting  equation  A>1  into  equation  A>10  gives 


*Adi2  ■  -  c‘u‘  ♦ 


(A.12) 


where 


6  > 


du 

7i  • 


Equation  A-12  is  die  staffing  pdnt  for  fimnulations  wUdi  employ  time>iiidependent  fandamefual 
aotmions.  The  treamient  of  dretmdent  term  on  die  left  hand  tide  of  dus  equation  diffois  in  die  schemes 
piDpoaed  in  die  literature  (Wtobd,  Teiles,  and  Bidibia  1986). 


S1EADY-STATE  SOLUnON  PROCEDURE 


To  flhistiaie  die  procedure  fiu' obtaining  a  solutim  to  equation  A>12  over  die  boundaiy,  consider  first 
die  simple  steady-state  case  in  wMdi  6  ■  0.  The  im^ral  term  on  die  left-hand  tide  of  equation  A-12 
becomes  zero  and  die  ftdlowing  equation  results: 


c4i‘  + 


(A-13) 
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Savoie  the  boundaiy  is  2-D,  uaooto,  nd  divided  into  N  dements  1^  s  series  of  stnught  line 
sefSBCnts  vdiose  ento  Ik  (» the  boundsry  carve,  tsid  siicoessivdy  join  one  snodier  to  f(wm  s  closed  dicuit 
dx»t  the  boundaiy.  The  points  at  the  middk  of  each  demem  ivhete  imknown  vdues  On  equatkn  A-13) 
are  considered  are  cdkd  "nodes."  Boundaiy  dements  of  diis  typt  are  called  "constant"  dements. 

In  this  case  of  constant  elements,  die  values  of  tt  and  q  are  assomed  to  be  constant  over  each  dement 
and  equal  to  die  vdue  at  die  mid-dement  node.  Eqaatkn  A-13  can  be  discretized  for  a  given  point  "i" 
befiue  ^iplying  any  boundaiy  oonditioas  as 

II  II 

+  ruq*dr-£  fr.qu"dr  (A-14) 

2  j.i  j!  j.j  J 

vriiere  the  pdnt  "i"  is  one  of  the  boundaiy  nodes  and  Tj  is  the  boundaiy  of  the  "j”  donem.  For  this  type 
of  element  0.e..  constant)  the  boundaiy  is  dways  smooth,  since  the  node  is  at  the  center  of  the  demeitt. 
Hence cl«  k. 

The  u  and  q  values  can  be  taken  out  of  die  int^rels  since  they  are  constam  over  eadidcmem.  They 
are  labelled  u*  and  ^  for  each  element  "j,"  so  etpudon  A-ld  becomes 

m 

jLuU  I  A^u^-  r  0«q^  .  (A-15) 

2  j.i  j-l 

vriiere 


ft**- J[q*dr  .  G««j[u*dr  .  (A-16) 

These  imegnOs  idak  the  "i"  node  vrirere  the  fundamental  solution  is  acting  to  any  other  node  "j." 

If  it  is  assumed  diat  die  podtion  i  can  also  vary  ftom  1  to  N,  coire^iondiiig  to  an  qifdication  of  the 
fundamental  sdution  d  eadi  node  successively,  dm  one  obtains  a  system  of  equations  resulting  &om 
applying  equation  A-IS  to  each  boundaiy  point 
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ftitnr.  If  COB  kis 


when  iwj 
when  i=j  * 

2 


dm  equatkn  A-16  can  be  written  as 

N  N 

E  H«u*  =  E  G«q^ . 
j»l  j.l 

This  set  of  equations  for  die  unknowns  and  q*  can  be  written  in  matrix  form  as 

HU  =GQ. 


(A-17) 


(A-lg) 


(A-19) 


where  H  and  are  two  N  x  N  matrices,  and  U,  U  tune  vectors  of  length  N. 

Snoe  there  are  values  of  u  and  N2  values  of  q  known  on  portions  Fj  and  r2  reflectively,  udiere 
N|  4  N2  «  N,  there  are  only  N  unknowns  in  die  system  of  equations  A>19.  To  introduce  these  known 
boundny  values  into  equation  A«I9,  the  system  must  be  rearranged  by  moving  columns  of  Hand  Qfrom 
one  side  to  die  odier.  When  all  die  unknowns  are  passed  to  the  left-hand  side,  the  system  can  be  written 

(A-20) 

iriiere  £  is  a  vector  of  unknown  (u's  and  q's)  boundary  values. 

Equation  A-20  can  be  solved,  and  tdl  the  boundary  values  wUl  dien  be  knowa  dus  infinmation, 
any  intenud  value  of  u  (n  its  derivatives  can  be  calculated  at  any  internal  poim  T  using  the  corresponding 
steady-state  fimn  of  equation  A-7, 


u‘=  f  qu*dr  -  f  uq*dr  .  (A-21) 

Jr  Jr 

Ifere  die  fundamental  solution  is  considered  to  be  acting  tm  an  internal  point  "i,"  and  all  the  values  of  u 
and  q  on  die  boundary  are  already  known. 

The  process  is  dien  one  of  numerical  iniegratioa  The  same  discretizadmi  is  used  for  die  boundaiy 
integrals 
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(A-22) 


N  N 

u‘-  I  Ciqi-  I  . 

j-1  j-1 

and  die  ooeffickma  and  fi**  must  be  calculated  again  fix  each  difierent  imemal  poim  L 

TIME  INDEPENDENT  FORMULATION 

To  proceed  with  die  fimnuladan  of  time-independent  fimdamental  atduttons  uaing  equation  A-12.  die 
ao-caUed  Dual  Reciprocity  BEM  developed  by  Wrobelndoo-wmkmsCWrobel  1988;  Wrobel,  Telles,  and 
BiebbU  1986)  wOl  be  outlined  here. 


The  function  d  at  any  pctint  inside  die  domain  O  is  ifiproximated  1^  a  set  of  N  (diosen)*  cooidinate 
functions  f*(x)  multiplied  by  ftmctioos  d^t).  ititich  are  unknowns. 


N 

d(jt.t)-  I  fKx)aKt). 

j-1 

Using  diis  iqjproximation.  die  domain  int^ral  on  the  left-hand  side  of  equation  A-12  becomes 

N 


(A.23) 


fu’ddQ-  £  dt^  ff^u* 
Ja  Ja 


da 


(A-24) 


It  is  furdier  assumed  diat  there  exists  a  function  i|H(x)  fix  each  function  f^x),  vdiere  the  two  functions 
are  related  by 


V*vJ(x)-fUx) 

Substituting  ecpiation  A-25  into  equation  A-24  gives 

N 


J^u'ddfl-  I  di^J^V*y^u* 


dO 


(A-25) 


(A-26) 


Introducing  eiqirestion  A-26  into  equation  A-12  yields  the  result 


*Cliooriin  dw  coowHimf.  finctiao  f‘  will  be  briefly  eddweeed  Uiei. 
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(A-27) 


N 


c*ii‘  + J^q’iidr- J^u*qdr«  £  idij|-J^tt*kV*i|rJdQ| . 

Now  applying  die  nanaftimnion  indicatBd  in  eqiution  A>10  to  the  domain  intend  gives  die  result 

cHi'+J^q’udr-jllu'qdr-  £  i|c‘y  ♦j|]q*\|r^dr-j^u*q^dr|&^  (A-28) 


^i^iere 


Ecpiaiion  A>28  involves  boundaiy  integrals  only. 

Wfobd  and  oo>woikei8  (Wiobel  1988;  Wiobel,  Telks,  and  Brebbia  1986)  found  it  convenient  to 
approximtte  die  variation  of  functions  u.  q.  y,  nd  q  witbin  each  bmmdary  element  using  a  unique 
(unqwdfied)  set  of  intetpolarion  functions  such  that 

u  -  O^u, 


q  ■ 


and 


(A-29) 


qi .  ©’’qi 

triiere  subset^  e  denotes  a  pafdcular  boundaiy  element 

Wiobel  (1988)  claims  dut  it  is  not  necessary  to  qjpiDximate  the  functions  y  and  q  dus  way  over  the 
boundaiy,  since  diey  are  known  once  the  set  of  functions  t  are  chosen.  However,  he  notes  dutt  this 
approxinuttion  dramatically  reduces  the  retpiired  boundary  integrations,  with  smne  sacrifice  in  accuracy. 
The  seriousness  of  this  sacrifice  is  tudmown,  but  should  be  exidoied. 
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l^xn  dividiof  the  bomdaiy  into  dements,  i^iplyiitg  the  discretized  versian  of  equstioo  A-28  to  each 
boundny  node,  and  taking  the  peeoeding  appioximaiioas  into  account,  the  Idlowing  system  of  cquatkms 
(witaen  in  matrix  fbim)  results: 


HU-OQ-I(HUv-Gil)a  ■ 

By  evaluating  the  discretized  vmsion  of  equation  A-23  at  all  boundary  nodes,  one  finds 

tf-Fft  . 


(A-3(» 


(A-31) 


Upon  inversion,  this  yields  die  result 


A  -  P‘  tl  . 


Substitudrig  eiqptession  A-32  into  equation  A>30  gives  the  matrix  form 

CtJ  +  aiJ  -GQ 


where 


C-  ^(Hv-Gn)P‘  . 


(A-32) 


<A-33) 


In  die  matrix  equation  A-33,  the  vectors  U  and  Q  each  have  the  dimension  equal  to  the  number  N  of 
boundary  nodes,  and  die  rrutrkes  C,  H  and  G  are  of  die  size  0>e.,  N  x  N)  of  standard  boundary  demod 
matrices. 


Wrobd  and  co-workers  (Wrobd  1988;  Wrobd.  Tdles,  and  Brebbia  1986)  furdier  utilized  a  simide, 
two-levd  time  integration  scheme  in  tbdr  work.  This  involved  using  a  linear  qiproximadon  for  U  and 
Q  vriddn  eadi  time  stq>  of  the  form 

U-(l-P)U-  +  pU-* 

Q-(1-P;Q"  +  PQ“*‘ 

(J  ■  J-  (ir*‘  -  U")  (A-34) 

At 

^’riiere  P  (O  <  P  <  1)  is  a  parameter  which  indicate  the  vdues  of  U  and  Q  bdween  time  levels  m  and 
m  1. 
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SulMlitutiiig  tliBte  approximations  iiao  equation  A-33  results  in  the  fidlowiqg 


right-hand  side  of  equation  A-35  is  known  at  time  m  At  By  introducing  the  boundaiy  conditions  at 
(m  -f  1)  ^  the  left^iand  side  of  equation  A-3S  can  be  reamnged,  and  die  system  of  equations  can 
>lved  for  this  time  step. 

SELECTION  OF  CXDORDINATE  FUNCTIONS 

T  selection  of  cooidinate  functions  f^x)  used  to  qjproximate  die  time  derivative  of  the  ^qxndett 
variaUe  u(x4:)  over  the  domain  O  is  qipaiendy  quite  impoitant  Ceitainly  it  needs  further  attention  in  any 
future  effort  to  apfdy  this  overall  BEM  to  die  barrel  heating  proUem. 

Wrobel  (193^11  conjectures  diat  die  choice  of  functions,  f^x).  should  conespond  to  die  behavior  of  die 
fundamental  solution  itrelf.  He  has  iqiparendy  had  some  success  in  2-  and  3-D  proUems  with  Cartesian 
coordinates  using  fonctiom  of  die  fonn 

fJ(x)  -  R(x.xJ)  .  (A-36) 

Here  R  rqiresents  die  Euclidean  distance  between  two  pmnts  in  space,  x*  are  coordinates  of  die  jdi  fixed 
point,  and  x  are  the  coordinates  of  any  point  in  the  domain  Q.  This  set  of  coordinate  functions  is  linearly 
indqiendent  as  Iraig  as  the  fixed  points  are  separate  (not  coincident).  The  corresprading  functims  y  and 
q  are 

yi(x)-iR»(x,xJ)  (A-37) 

a 


iju* 


C-(l-p)H  U*  +  (1-P)GQ- 


(A-35) 


and 


qJ(x)-i!iR’(x,xJ) 

a 


dR 

■ST’ 


(A-38) 
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a  hai  die  value  of  9  or  12  for  2-  or  3-D  probtems,  le^Kcdvdy. 


For  axisymuietric  piDbleina,  Wrobel  TeUes,  and  Brebbia  (1986)  have  suggested  die  use  of  oooidiiiatB 
funcdoos  of  the  fimn 


fJ(x)-R'(x.xJ)Il-^J  , 


(A-39) 


v^ieie  (or  r)  is  die  distmce  horn  a  fixed  pdnt  (or  soy  point  in  die  domain  Q)  to  the  axis  of  levtdudon, 
and  is  die  Euclidean  distance  between  points  defined  in  the  (r^)  fdane  Gx..  on  the  generating  pine 
0  «  constant).  In  this  case,  the  eiqnessions  for  y  nd  are  of  the  fimn 


y(x)-.^R'*(x.xJ)  (A-40) 

and 

tii(x)-^R'‘(x.x^)^.  (A-41) 


Wrobd  (1988)  also  notes  the  following  items: 

(1)  The  functions  F  should  be  lineady  independent  for  die  matrix  F  to  have  an  inveise. 

(2)  Once  die  family  F  is  chosn,  it  is  usually  easy  to  determine  die  fonctirais  and  iV. 

(3)  The  position  and  number  of  fixed  points  is  ctmveniendy  chosen  to  coincide  with  the  boundaiy 
nodes. 

(4)  The  inclusion  of  an  internal  poim  associated  with  f  =  amstam  may  be  required  to  better  simulitte 
the  heating  up  of  the  entire  body  by  a  cemstant  value. 
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Time  weie  two  miyor  enois  in  die  originai  XBR-2D  ideaae  dm  weie  discovered  and  pointed  out  to 
Vetitay  by  Anny  Research  Laboratory  (ARL)  pefanmel  near  the  end  of  this  effoii  Hguie  B-1  Ohistnies 
the  results  fram  several  nms  with  die  original  code.  These  curves  rqaeseitt  die  surface  temperature  rise 
at  die  42’in  location  of  a  155-inffi  howitzer.  The  only  unidicit  input  differences  between  the  seven  nns 
are  die  locations  of  the  left  and  right  hand  computational  boundaries.  The  resulting  temperatures  are  very 
different  from  each  other.  It  was  this  finding  that  called  atiendtm  to  the  errors  in  coding. 


Hguie  B-1.  Orieinal  XBR-2D  runs  with  various  computational  regions. 


One  of  die  inteipolatitm  functions  in  the  code,  whidi  converts  the  geometric  parameters  fiom  die  input 
file’s  axial  positions  to  that  of  the  individual  grid  locations,  had  a  mistake  in  die  foimidatioiL  The  left 
side  boundary  was  always  assigned  the  geometric  values  of  the  first  input  value  (inside  and  outside 
diameters),  no  matter  where  this  boundary  was  located.  If  the  left  side  boundary  was  not  equal  to  die  first 
geometric  location  in  the  irqxit  file,  the  wrong  geometry  was  used  for  that  point.  This  has  been  corrected 
in  the  new  version  of  die  code. 

A  more  significant  error  present  in  the  original  code  was  the  treatment  of  the  diermal  boundary  layer. 
The  diickness  of  the  boundary  layer  was  always  set  to  zero  at  the  left  hand  boundary.  This  is  not  correa 
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fiir  a  kfk  ikte  teundary  that  does  not  suit  at  the  rear  of  the  gun.  lUs  enor  resulu  in  a  much  tWnner 
boundaiy  layer  than  is  ooirect  if  the  left  side  boundary  is  not  at  the  rear  of  the  gun. 


Ihis  enor  was  corrected  by  adding  subfoutine  FZE^IN  t^iich  calculates  die  thickness  of  the  diennal 
boundary  layer  up  to  the  left  side  boundary  and  uses  diis  as  an  initial  conditioo  for  the  calculations.  The 
test  cases  nm  in  Figure  B-1  were  rqieated  aftm  die  correcdons  were  made  and  are  piesenied  in 
Hgure  B-2. 


The  agreement  between  die  curves  is  reasonable  considering  the  wide  range  of  qntial  resoludons  diat 
are  iqxesented  in  die  figure.  The  first  curve  defined  the  simulation’s  100  grids  to  represent  the  emiie 
band,  dms  ftndng  die  axial  grid  qiadiig  (Az)  to  be  3.13  ia  Oonverady,  die  last  amdadon  presented 
defines  the  code’s  100  grids  over  a  mete  10  in.  making  the  axial  qwcing  0.1  in.  For  simulations 
conducted  Hcmg  previous  verdons  of  die  XBR-2D  code,  dre  results  will  change  slisftfiy  if  the  left 
cooqxitational  boundary  was  not  set  to  zero.  Simuladons  widi  a  left  boundaiy  at  z  s  o  should  be 
essentially  unchanged  finm  the  previous  version  of  the  code. 
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APPENDIX  C: 

FURTHER  IMPROVEMENTS  TO  THE  XBR-2D  CODE 


iNnminMLLY  LEFT  BLANK. 
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The  OiliiDil  XBR-2D  omIb  cnwained  tone  am  Id  cocttig  nl  in  dieoiy.  llieae  oion  wea 
ooimed  «Ml  iaploMaied  in  voiion  Tte  addition  of  muUiite  grids  in  die  chmelaya,  if  ttnc 
liga  beoomi  tuck,  was  also  added  as  a  new  feaniie  of  venion  1 JX). 

FINITE  OnERENCZ  EQUATION 

Tlie  original  riniie  dUfennoe  mnation  in  XBR  was  oonect  if  die  geometric  factm  was  oonatanr 
dumi^toitttliettdGfaieas.  However,  the  code,  in  fitting  die  gridding  to  die  band  geometty,  modified  the 
grids  andi  that  the  geometric  fidor  between  the  flat  and  second  grid  was  diffinent  fram  the  geometric 
ftcior  between  the  rest  the  grids. 

To  coned  dds  error  and  to  accommodate  the  modification  which  allows  for  muldple  grids  throng 
die  chrome  layer,  the  finite  difference  equation  was  modified  to  allow  for  a  general  gridding  sdieme  with 
different  radial  grid  yicing  in  dther  dtaecdoa  The  Ar  is  computed  for  each  radial  position  and  it  is  this 
value,  not  the  geometric  fector.  wUdi  is  used  in  the  finite  difference  equation. 

For  die  Fourier  eolation  of  heat  oanduction  in  c^indrical  coordinates. 


«»r  1  «T  1  «  /A  8T.  1  5T 

8?  7  Tr  -X:  ‘51  "  o  TT  ‘ 

Hie  generalized  finite  difference  equation  fix'  variatde  Ar  is: 


(C-1) 


a  At 


T  -?T  -*-T  -  IT  -T 


i(4r..4rj 


[Tg..-Ty.J.^[T^,j-2Ty*T.,.J 


(C-2) 
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The  fUie  dilfeicaoe  equttioo  ibr  varieble  /h*  at  the  duDoe-aied  inteifMe  is: 


(Ar**+Ar*.)  L 


-At  ) 

®.Tu*j  -  2«,«Ty  ♦  «,Tg-j  -  Pu**  "  ^U-J 


*  tuT,  •  (C-3) 

These  two  eqiiatiaiis  weie  mqdemented  in  XBR  venim  1J20  to  replace  flie  original  finite  difference 
equations. 


RADIAL  GRIDDING 


In  the  (Miginal  venions  of  XBR-2D.  the  number  of  radial  grids  was  computed  based  on  die  local 
diidcness.  This  could  result  in.  for  certain  geometries,  different  numbers  of  radial  grids  as  the  thidmess 
varied  axially.  In  such  a  situation,  the  axial  heat  conductton  calculations  were  invalid.  This  was  oonected 
by  computing  die  number  of  grids  required  for  the  diinnest  part  of  die  barrel  and  stretdung  that  number 
of  grids  to  fit  the  thickness  at  odier  axial  posidons. 

Two  odier  ahematives  were  tested:  (1)  basing  die  number  of  grids  tm  the  thidcest  part  of  die  barrel 
and  shrinking  the  grids  to  fit  into  the  thinner  parts,  and  (2)  using  an  average  thickness  and  both  shrinking 
and  stretdung.  vdiere  apprc^ate,  to  fit  die  barrel  geonretry.  Both  of  these  methods  caused  the  grid  sizes 
to  become  extremdy  small  in  the  diinnest  secdmi  of  die  barrel  when  grid  dirinking  was  restdved.  This 
resulted  in  a  drasdc  reduction  in  the  dme  step  and  a  corresptxidingly  large  increase  in  execution  time. 

mulufle  grids  in  the  chrome  layer 

XBR-2D  originally  fixed  the  th«ckne«?  of  die  first  grid  to  the  thickness  of  the  duome  layn,  if  sudi 
a  layer  was  pre*”**  If  die  duome  layer  becomes  thick,  or  XBR  is  used  to  modd  other  types  of  liners. 
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tbe  dHanetoyer  mult  be  modeled  widi  more  than  one  grid.  XBR  was  modified  ao  that,  if  a  chiome  layer 
is  inesent.  die  fim  grid's  dddmess  is  based  on  die  rime  step.  If  the  chrome  thickness  is  greater  than  diis 
dddmess.  dien  the  ddcknem  of  the  second  grid  is  computed  based  on  the  geometric  faatot  and  this  prooem 
is  continued  until  die  chrome  thickness  is  readied.  The  exact  thickness  of  the  dirome  is  then 
acoompUshed  by  shrinldiig  the  grids  to  this  thickness.  The  inteifitte  between  the  dirome  and  sted  » 
handled  by  Equation  2. 

THERMAL  BOUNDARY  LAYER 

The  dwimal  boundary  layer  is  calculated  using  the  method  described  by  Stratford  and  Beavers  (S&B) 
in  1961.  In  the  translation  fiom  the  prqier  to  the  code  several  mistakes  were  made.  The  code  originally 
had  die  Mach  number  wdghting  fiinctkm,  P,.  defined  as: 

P,-M^  1^1  (y  - 1)  M*j^  .  (C-4) 

After  reviewing  die  S&B  method,  it  was  determined  diat  the  equation  for  the  Madi  number  wei^iting 
fimctitm  should  be: 


M 


>r-« 

ray=n 


l  +  ^(y-l)M* 


.us 


(C-5) 


The  term  r*^  was  S&B's  correction  for  axisymnutric  flow.  The  other  corrections  involved  moving 
the  Mach  number  inside  the  brackets  and  correcting  and  moving  the  second  term  to  the  (tenominator  of 
the  exptessitML  This  corrected  equation  now  agrees  with  the  method  described  by  S&B  in  their  1961 
paper. 

MODIFICATIONS  TO  THE  RESTART  OPTION 


The  miginal  coding  of  the  restart  option  required  the  simulation,  when  restarted,  be  run  with  the  same 
cmifiguration  and  timing  as  the  original  run.  This  was  slightly  modified  in  the  new  version  of  the  code. 
The  code  now  writes  the  time  into  the  file  LASTTEMP.OUT.  When  the  restart  option  is  diosen,  the 
simulation  uses  this  time  as  the  starting  point  This  allows  simulations  of  various  firing  rates  to  be  rutt 
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NEW  INPUT  FE£  FORMAT 


Two  mote  ou^  flags  woe  added  to  die  iquit  flk  for  graphical  output  files,  and  an  additional  line  was 
added  where  ttie  barrel  and  chrome  material  properties  must  be  input  Ftdlowing  is  a  snqjle  irqwt  file 


r  S«Npl«  xm-2D  Input  fll*  for  •  2S>«b  Inrrol  (7S0  Shota/aln) 


t  tablont  Tamp  (R>,  Pr  (pal),  gaa  «•!  (In/a) 

;  :  . :  : 

;  >  . 

;  I  :  : 


S30.0  14.70  0.00 


i  tattol  Coaaatry  < 


Bchrcn  (In).  tiMxv  (na),  aeon,  aehaa  (In),  nax  (In) 
t  :  s  : 

. :  :  : 

. . :  : 

:  I 


t  nloe  .  catart  (In),  cond  (In) 


;  : 

:  : 

t  i 

:  : 

;  : 

>  : 

i  t 

a  a  «  a 

at  : 

;  : 

s  : 

:  : 

t 

0.005 

4080.0  0.100 

4.60  73.5 

20 

0.0 

10.0 

ndtz. 

cdb  (lindtz) 

—  no  of  axial 

loeatlona,  ndtz  loca 

(■ax 

-  8) 

1 

2  3 

4  5 

6 

7 

8 

7  0.0 

4.0  13.0 

23.0  43.0 

«5.0 

75.0 

dbar 

(l:ndtz) 

—  barzol 

o.d. 

at  ndtz 

loeatlona  (1 

1 

2  3 

4  5 

6 

7 

8 

2.75 

2.75  2.25 

2.00  1.90 

1.75 

1.65 

dbor  (l:ndtx)  —  barral  l.d.  at  ndtt  loeatlona  (in) 

;12345678 
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1.0>4  l.OM  1.024  1.024  1.024  1.024  1.024 
I 

}  OMixad  0«tp«t  lOMtlOMt 

I.— ..II.  I  -  - I  . . . . 

I  Bkts.  Mat  (liBhts)~  M  of  leeotioaa  for  output,  ohta  loeotloua  (la) 
I  (MB 
t 

;1234Sf7t 

t 

•  4. 425  4.425  4.425  4.425  7.0  7.0  7.0  7.0 


rout  (liahtt)  —  radial  dapth  at  ahti  loeatloM  (la) 
12345474 

0.0  0.005  0.025  0.5  0.0  0.005  0.02  .5 


; 

;  Flrlag  Data: 

; - - - -  ■■ 

;  arounda,  n)countl,  nkountZ,  daitv  (a),  laataliot 


;  . .  t  I  t 

t  I  . .  s  I 

>  Si  . s  : 

)  1  I  s  . 

;  t  I  s  : 


1  10  50  i.a>-s  1 

; 

l  Ooalrod  Output  Units  and  riloss 

; - — - - -  ,  . 

;  lunlt  (0  for  K,  1  for  dog  R) 

;  ...  ItOBp  (0  for  at  ealeulatod,  1  for  tsap  rise), 

;  s  lout  (1:9)  - 

;  :  .  scroon  output 


;  :  : 

;  :  :  . tlnglo  allot  prolM  ttaporaturo, 

;  I  :  :  . tlnglo  t)iot  Intldo  turfaco  taaporatoro, 

;  :  :  :  :  . tlnglo  tlwt  outtldo  turfaco  ttaporaturo, 

;  :  :  :  :  :  . burst  flro  probe  toaporaturo, 

;  :  :  :  i  :  :  .  burst  flra  Intldo  turfaco  toaporaturo, 

.  burst  flro  outtldo  turfaco  toaporaturo. 


;  :  I  :  :  i  :  :  :  . )ioat  trantfor  eoofflclont, 

.  offoetlvo  gat  toaporaturo  history 

;  i  :  s  I  :  :  :  1  s  1 

0002020100000 

J  •  • 

t  : . grap)ilcal  output  fllo 
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•ural 

k 


graphleal  output  fllo 


■orrol  ChroM  ChroM 

a  k  a 


4.C4t2  O.OISC  ».3S$i  0.03131 


Surface  Temperature  (K) 


rface  Temperature  (K) 


lUe^VtBiiiXwT 


Te«Na4:  Examine  the  Mach  Dumber  profile  down  the  barrel 


Figure  C-4.  Madi  number  profiles  vs  time  and  position  for  tfie  155-mm  howitzer  test  case. 

I 

I 
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TottNo.  S:  Apply  a  oonnat  liB«  flux  to  the  bon  nifitte  nd 


ttie  leiuits  to  an  aulyikal 


Figure  C-S,  Comparison  of  XBR  witti  a  constant  heat  flux  to  an  analytical  solution. 


A  heat  flux  of  IJOOOjOOO  watts  was  qiplied  to  the  suiftce.  The  analytical  solution  was  obtained  from 
tibe  equation: 


No.  or 


Nouor 


fiBBli  fiiwilnB 
2  AAmMmwn 

TDChoiCOl  C'minr 

ATIN:  OnC-DDA 
SttdoB 

AkxMMfeia.VA  22304-6145 

1  OmimiMW 

U^.  Anny  Mtterid  Cnmnwid 
ATIN:  AMCAM 
5001  EiMohower  Ave. 

Alexandria.  VA  22333-0001 

1  DinciQr 

U.S.  Anny  Reaeareli  Labomory 
ATIN:  AMSRL-OP-a-AD. 

Tech  Publirtag 
2800  Powder  MOlRd. 

AdeIphi.MD  20783-1145 

1  Dhecior 

U.S.  Anny  Reaewcfa  Labocaiaty 
ATIN:  AMSRLOP-a-AD. 

Reoanb  Management 
2800  Powder  MiURd. 

AdeIphi.MD  20783-1145 

2  Comnunder 

U.S.  Anny  Annamem  Reaevch. 

Devdoproetu.  and  Engineering  Center 
ATIN:  SMCAR-TDC 
Picatinny  Anenal.  NI  07806-5000 

1  Direcaor 

Laixmoiy 

U.S.  Anny  Amument  Raearch. 

Development,  and  Engineering  Cenia 
ATIN:  SMCAR-CCB-TL 
Waiervliet.  NY  12189-4050 

1  Diiecior 

U.S.  ^my  Advanced  Systems  Research 
and  Ani^fsis  Office  (ATCOM) 

ATIN:  AMSAT-R-NR.  M/S  219-1 
Ames  Research  Center 
Moffett  Field.  CA  94035-1000 


1  Connnander 

U.S.  Amy  Muaile  Command 
ATIN:  AMSMI-RIMS-R  (DOC) 

Redstone  Aiaend.  AL  35898-5010 

1  Commander 

U.S.  Amy  Tank-Antomoiive  ComnMnd 
ATIN:  AMSTA-JSK  (Armor  Eng.  Br.) 
Wanen.MI  48397-5000 

1  Diiecior 

U.S.  Amy  TRADOC  Analysis  Command 
ATIN:  ATRC-WSR 

White  Sands  Miasik  Rai«e.  NM  880Q^55Q2 

(□■■.■air)!  Commandant 

U.S.  Amy  bfuMiy  School 
ATIN:  ATSH-CD  (Secnrity  Mp.) 

Fon  Benning.  GA  31905-5660 

(u*^*h)l  Commandant 

U.S.  Amy  Infimiry  School 
ATIN:  ATSH-WCB-O 
Fort  Benning.  OA  31905-5000 

1  WL/MNOI 
EghnAFB.FL  32542-5000 

Aberdeen  Provine  Oround 

2  Dir.USAMSAA 
ATIN:  AMXSY-D 

AMXSY-MP.  H.  Coben 

1  Cdr.  USATECOM 
ATIN:  AMSTE-TC 

1  Dir.  USAERDEC 
ATTN:  SCBRD-RT 

1  Cdr.  USACBDCOM 
ATTN:  AMSCB-Cn 

I  Dir.USARL 

ATTN:  AMSRL-SL-I 


5  Dir.  USARL 

ATTN:  AMSRL-OP-AP-L 


51 


No.  oi 

Qsm  OrgiBiaUQP 

1  Chainnaii 

DOD  Eiqrfosives  Safety  Board 
Romn  856-C 
Hofiraao  Bldg.  1 
2461  Eiseahower  Avenue 
Aleundria.VA  22331-0600 

1  Headquarters 

U.S.  Army  Materiel  Command 
ATTN:  AMCICP-AD,  M.  Fisetle 
SOOl  Eisenhower  Ave. 

Alexandria,  VA  22333-0001 

1  U.S.  Army  Ballistic  Missile 

Defense  Systems  Command 
Advanced  Technology  Center 
P.O.  Box  ISOO 
Huntsville.  AL  3S807-3801 

1  Department  of  die  Army 

Office  of  the  Product  Manager 
155mm  Howitzer.  M109A6,  Paladin 
ATIN:  SFAE-AR-HIP-IP, 

Mr.  R.  De  Kleine 
Picatinny  Arsenal,  NJ  07806-5000 

3  Project  Manager 

Advanced  Field  Artillay  System 
ATTN:  SFAE-ASM-AF-E, 

LTC  A.  EUis 
T.  Kuriata 
J.  Shields 

Picatinny  Arsenal,  NJ  07801-5000 

1  Project  Manager 

Advanced  Field  Artilleiy  Syston 
ATTN:  SFAE-ASM-AF-Q,  W.  Warren 
Picatinny  Arsenal,  NJ  07801-5000 

2  Commander 

Production  Base  Modonization  Agency 
U.S.  Army  Armament  Research, 
Development,  and  Engineering  Center 
ATTN:  AMSMC-PBM,  A.  Siklosi 

AMSMC-PBM-E,  L.  Laibson 
Picatinny  Arsenal,  NJ  07806-5000 


No.  of 

Copies  nrg«ni7^rinn 

4  rcO-Armaments 

Project  Manager 
Tank  Main  Armament  System 
ATTN:  AMCRM-TMA 

AMCH4-TMA-105 
AMCnd-TMA-120 
AMCPM-TMA-AS,  H.  Yuen 
Picatinny  Arsenal,  NJ  07806-5000 

4  Commander 

U.S.  Army  Armament  Research, 
Development,  and  Engineering  Center 
ATTN:  SMCAR-CCH-V, 

C.  Mandala 
E.  Fennell 

SMCAR-CCH-T,  L.  Rosendorf 
SMCAR-CCS 

Picatinny  Arsenal,  NJ  07806-5000 

18  Commander 

U.S.  Army  Armament  Researeh, 
Development,  and  Engineering  Center 
ATTN:  SMCAR-AEE,  J.  Lannon 
SMCAR-AEE-B, 

A.  Beardell 
D.  Downs 
S.  Einstein 
S.  Westley 
S.  Bernstein 
J.  Rutkowski 

B.  Brodman 
P.  O'Reilly 
R.  Cirincione 
P.  Hui 

J.  OReUly 
SMCAR-AEE-WW, 

M.  Mezger 
J.  Pinto 
D.  Wiegand 
P.  Lu 

C.  Hu 

SMCAR-AES,  S.  Kaplowitz 
Picatinny  Arsenal.  NJ  07806-5000 

I  Commander 

U.S.  Army  Armament  Research, 
Development  and  Engineoing  Center 
ATTN:  SMCAR-HFM,  E.  Barrieres 
Picatinny  Arsenal,  NJ  07806-5000 
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No.  of 

QmnixitigD 

8  Commander 

U.S.  Army  Armament  Research, 
Development  and  Engmeering  Center 
ATTN:  SMCAR-FSA-F,  LTC  R.  Riddle 
SMCAR-FSC,  G.  Fddinand 
SMCAR-FS.  T.  Gora 
SMCAR-FS-DH.  J.  Feneck 
SMCAR-FSS-A. 

R.  Kopnunn 
B.  Machek 
L.  Finder 

SMCAR-FSN-N.  K.  Chung 
Picatinny  Arsenal,  NJ  07806-5000 

3  Director 

Benet  Weapons  Laboratories 
ATIN:  SMCAR-CCB-RA, 

G.P.  O'Hara 
G.A.  Pflegl 

SMCAR-CCB-S,  F.  Reiser 
Watervliet,  NY  12189-4050 

2  Commando' 

U.S.  Army  Research  Office 
ATTN:  Technical  Library 
D.  Mann 
P.O.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 

1  Commando,  USACECOM 

R&D  Technical  Library 
ATTN:  ASQNC-ELC-K-L-R, 

Myo  Cento 

Fort  Monmouth,  NJ  07703-5301 

1  Commandant 

U.S.  Army  Aviation  School 
ATTN:  Aviation  Agency 
Fort  Rucko,  AL  36360 

1  Program  Manago 

U.S.  Tank-Automotive  Command 
ATTN:  AMCPM-ABMS,  T.  Dean 
Warren,  MI  48092-2498 


No.  of 

Copies  nry«iiiMtion 

1  Project  lago 

U.S.  TanL -.'Automotive  Command 
Fighting  Vehicle  Systems 
ATTN:  SFAE-ASM-BV 
Warren,  MI  48397-5000 

1  Project  Manager,  Abrams  Tank  System 

ATTN:  SFAE-ASM-AB 
Warren,  MI  48397-5000 

1  Director 

HQ,  TRAC  RPD 

ATTN:  ATCD-MA 

Fort  Monroe,  VA  23651-5143 

1  Commander 

U.S.  Army  Belvoir  Research  and 
Development  Center 
ATTN:  STRBE-WC 
Fort  Belvoir,  VA  22060-5006 

1  Director 

U.S.  Army  TRAC-Ft  Lee 
ATTN:  ATRC-L,  Mr.  Cameron 
Fort  Lee,  VA  23801-6140 

1  Commandant 

U.S.  Army  Command  and  General 
Staff  College 

Fort  Leavenworth,  KS  66027 

1  Commandant 

U.S.  Army  Special  Warfare  School 
ATTN:  Rev  and  Tmg  Lit  Div 
Fort  Bragg,  NC  28307 

1  Commander 

Radford  Army  Ammunition  Plant 
ATTN:  SMCAR-QA/HI  LIB 
Radford,  VA  24141-0298 

1  Commander 

U.S.  Army  Foreign  Science  and 
Technology  Center 
ATTN:  AMXST-MC-3 
220  Seventh  Street,  NE 
Charlottesville,  VA  22901-53% 
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No.  of 

CflBkft  OfitiT-Mtiftn 

2  ComimndaBt 

U.S.  Anny  Keid  AitiUery  Center  and 
School 

ATTN:  ATSF-CO-MW.  E.  Dublisky 
ATSF-CN,  P.  Gross 
R  Sill.  OK  73503-3600 

1  Commandant 

U.S.  Army  Armor  School 

ATTN:  ATZK-CD-MS,  M.  Falkovitch 

Armor  Agency 

FwtKnox,  KY  40121-5215 

1  U.S.  Army  Research  Development  and 

Standardization  Group  (UK) 

PSC  802  Box  15,  Dr.  Roy  E.  Richenbach 
FPO  AE09499-1500 

2  Conunander 

Naval  Sea  Systems  Command 
ATIN:  SEA62R 
SEA  64 

Washington.  DC  20362-5101 

1  Commander 

Naval  Air  Systems  Command 
ATTN:  AIR-954-Tech  Library 
Washington.  DC  20360 

4  Commander 

Naval  Research  Laboratory 
ATTN:  Technical  Library 
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2  Commander 

Naval  Surface  Warfare  Center 
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Naval  Surface  Warfare  Center 
ATTN:  T.C.  Smith 
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STE9 

EglinAFB.FL  32542-5910 

1  WL/MNSH 
ATTN:  R.  Drabczuk 
EglinAFB.FL  32542-5434 

2  NASA  Langley  Research  Center 
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1  Rensselaer  Polytechnic  Institute 
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Davidson  Laboratory 
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Aerospace  Engineering 
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1  University  of  Utah 
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1  AFELM.  The  Rand  Corporation 

ATTN:  Library  D 
1700  Main  Street 
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1  Arrow  Technology  Associates,  Inc. 

ATTN:  W.  Hathaway 

P.O.  Box  4218 
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3  AAI  Corporation 
ATTN:  J.  Hebert 

J.  Frankie 
D.  Cleveland 
P.O.  Box  126 
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2  Alliant  Techsystems,  Inc. 

ATTN:  R.E.  Tompkins 

J.  Kennedy 
7225  Northland  Dr. 

Brooklyn  Park,  MN  55428 

1  AVCO  Everett  Research  Laboratory 

ATTN:  D.  Stickler 
2385  Revere  Beach  Parkway 
Everett,  MA  02149-5936 

1  Genera]  Applied  Sciences  Lab 

ATTN:  J.  Erdos 
77  Raynor  Ave. 

Ronkonkama,  NY  1 1779-6649 

1  Martin  Marietta 

Tactical  System  Department 
ATTN:  J.  Mandzy 
100  Plastics  Ave. 
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D.A.  Worrell 
WJ.  Worrell 
C.  Chandler 

Radford,  VA  24141-0299 


No.  of 

OrgMiizarion 
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1  Hercules,  Inc. 
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ATTN:  R.  Cartwright 
100  Howard  Blvd. 
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ATTN:  B.M.  Riggleman 
Wilmington,  DE  19894 

1  MBR  Research  Inc. 
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601  Ewing  St.,  Suite  C-22 
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ATTN:  F.E.  Wolf 
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1048  South  St. 
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3  Rockwell  Intnnational 
Rocketdyne  Division 
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J.  Flanagan 
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2  Rockwell  International  Science  Center 
ATTN:  Dr.  S.  Chakravarthy 

Dr.  S.  Palaniswamy 
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ATTN:  M.  Palmer 
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1  Southwest  Research  Institute 

ATTN:  J.P.  Riegel 
6220  Culebra  Road 
P.O.  Drawer  28S10 
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1  Sverdrup  Technology,  Inc. 

ATTN:  Dr.  John  Dcur 
2001  Aerospace  Parkway 
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ATTN:  R.  Wilier 
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Elkton.  MD  21921-0241 
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1  Universal  Propulsion  Company 

ATTN:  HJ.  McSpadden 
25401  North  Central  Ave. 

Phoenix,  AZ  85027-7837 

1  SRI  International 
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333  Ravenwood  Avenue 
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1  Cdr,  USACSTA 

ATTN:  STECS-PO/R.  Hendricksen 
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1  Erast-Mach-Institut 
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Germany 
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England 
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2  Institut  Saint  Louis 

ATTN:  Dr.  Marc  Giraud 
Dr.  Gunther  Sheets 
Postfach  1260 
78S8  Weail  am  Rhein  1 
Germany 

1  Explosive  Ordnance  Division 

AT^:  A.  Wildegger-Gaissmaier 
Defence  Science  and  Technology 
Organisation 
P.O.  Box  1750 

Salisbury,  South  Australia  5108 

1  Armaments  Division 

ATTN:  Dr.  J.  Lavigne 
Defence  Research  Establishment  Valcartia- 
2459,  Pie  XI  Blvd.,  North 
P.O.  Box  8800 

Courcelette,  Quebec  GOA  IRO 
Caiuula 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  puUishes.  Your 
comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  ARL  Report  Number  A»L-CR-126 _ Date  of  Report  April  1994 _ 

2.  Date  Report  Received _ 

3.  Does  this  report  satisfy  a  need?  (Grmmem  on  purpose,  related  project,  or  odier  area  of  interest  fw 

which  the  report  will  be  used.) _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved, 
operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  dianged  to  improve  future  reports?  (Indicate 
changes  to  organization,  technical  content,  format,  etc.) _ _ 


Organization 

CURRENT  Name 

ADDRESS  _ 

Street  or  P.O.  Bex  No. 

City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address 
above  and  the  Old  or  Incorrect  address  below. 


Organization 

OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 

City,  State,  Zip  Code 

(Remove  this  sheet,  fold  as  indicated,  tape  closed,  and  mail.) 
(DO  NOT  STAPLE) 


